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Abstract 

For slowly varying fields the vacuum functional of a quantum field theory may be 
expanded in terms of local functionals. This expansion satisfies its own form of the 
Schrodinger equation from which the expansion coefficents can be found. For scalar 
field theory in 1+1 dimensions we show that this approach correctly reproduces 
the short-distance properties as contained in the counter-terms. We also describe 
an approximate simplification that occurs for the Sine-Gordon and Sinh-Gordon 
vacuum functionals. 



1 Introduction 

Whilst asymptotic freedom has lead to an accurate determination of the Lagrangian of the 
Standard Model of particle physics from high energy experiments there are few analytical 
tools enabling us to calculate with that Lagrangian at low energies where the semi-classical 
approximation is no longer valid. For example, the eigenvalue problem for the Hamilto- 
nian of Yang-Mills theory cannot be solved in a semi-classical expansion because the 
renormalisation group implies that the energy eigenvalues depend non-perturbatively on 
the coupling. Consequently the computation of the hadron spectrum can only be done 
numerically. In ordinary quantum mechanics there are many ways to tackle this problem 
which are not widely used in a field theory context, but which, if suitably generalised 
might allow non-perturbative methods to be developed for field theory. The oldest of 
these is the Schrodinger representation, (see Q-|2l| for applications to field theory, and 



2"3R-|p9fl for applications to the Wheeler-de Witt equation). 

In the Schrodinger representation the vacuum, \E ), of a scalar quantum field theory 
is represented by the functional ((p\Eo) = expW[<p], where (tp\ is an eigenbra of the field 
operator 0(x) at fixed time, belonging to eigenvalue <£>(x). In general W is non-local, but 
if y?(x) varies slowly on the scale of the inverse of the mass of the lightest particle, mp 1 , 



it can be expanded in terms of local functionals, [p0 |, for example in 1+1 dimensions 



w = j dxJ2B 10 .,MxY°^(x) 



h ..^ n) {x) jn 



The coefficents Bj __j n are constant, assuming translation invariance, and finite as the ultra- 
violet cut-off is removed, ||. Particle structure is characterised by length scales smaller 
than tuq , so this simplification in W[(p] does not appear useful, however, a knowlege of 



this local expansion is sufficent to reconstruct W for arbitrary <^?@]- This is because if 
W[(p] is evaluated for a scaled field, <p s (x) = (p(x/y/s), it extends to an analytic function 
of s with cuts restricted to the negative real axis, so that Cauchy's theorem can be used 
to relate the large-s behaviour (when (p s is slowly varying) to the s — 1 value: 

W[<p] = lim — / -^-e^ s -^W[(p s }. (2) 

LhPJ a^oo 2ni J\s\=oo s - 1 [ ^ J V ; 

The exponential term removes the contribution of the cut as A — > oo. The vacuum 
functional satisfies the Schrodinger equation from which the coefficents Bj __j n can be found 
in principle, however care must be exercised because this equation depends explicitly 
on short-distance effects via the cut-off, whereas the local expansion is only valid for 
fields characterised by large length-scales, so we cannot simply substitute ([!]) into the 
Schrodinger equation and expect to be able to satisfactorily take the limit in which the 
cut-off is removed. However, we can again exploit Cauchy's theorem to construct a version 
of the Schrodinger equation that acts directly on the local expansion by considering the 



effect of a scale-transformation on the cut-off, as well as on the field, [30]. The Hamiltonian 
with a scaled cut-off acting on the vacuum functional evaluated for the scaled field again 
extends to an analytic function with cuts on the negative real axis. This enables the 
limit in which the short- distance cut-off is taken to zero to be expressed in terms of large- 
distance behaviour described by the local expansion for W[y]. This leads to an infinite 



set of algebraic equations for the coefncents Bj „j n . By truncating the expansion the 
Schrodinger equation offers the possibility of solution beyond perturbation theory in the 
couplings, however, before this is attempted it is essential to show that this formulation is 
capable of reproducing the results that can be obtained using the standard approach of the 
semi-classical expansion and Feynman diagram perturbation theory. In particular, since 
the method consists of building states out of their large distance properties it is important 
to show that it gets right the short- distance behaviour as contained in the counter-terms 
of the Hamiltonian. The purpose of this paper is to demonstrate that this short-distance 
behaviour is correctly reproduced by our approach to the Schrodinger equation in which 
we build the vacuum state from its large- distance behaviour, and that, at least to low 
orders, the resulting local expansion coincides with the Feynman diagram calculation of 
the vacuum functional. 

2 Semi-classical analysis of local expansion 

The classical Hamiltonian of ip 4 theory is J dx(^ (n 2 + <p' 2 + m 2 ip 2 ) + jg<p 4 ) in 1+1 di- 
mensions. In the Schrodinger representation the canonical momentum is represented by 
functional differentiation ft = —ihS/Stp(x), so the kinetic term leads to the product of two 
functional derivatives at the same point which we regulate by introducing a momentum 
cut-off p 2 < 1/e. The couplings must consequently be renormalised so that the Hamilto- 
nian has a finite action on the vacuum. Writing the vacuum functional as exp (W[<£>]/7i) 
gives the Schrodinger equation as lim e j F e [ip] = where 

Tl AW I ,,J l i ( SW Y + ^2 + M 2 (( .y | + 9,4 _ ,:,,., | ,;:,, 



\*- W + J <*(*(-{ 



and 



2 J I 2 I \5<p J w / 4! 

A e = fdxdy f JP e »(*iO - 6 [ = f dp2ir * (4) 

J J P 2 <i/e in dip{x)d<p{y) J P 2 <i/e d(p{—p)dip{p) 

where <p(p) = f dx<p(x) exp(— ipx). In perturbation theory the only divergent diagrams 
with external legs are tadpoles, and these can be removed by normal ordering the Hamil- 
tonian. This enables the e-dependence of the parameters to be calculated exactly as 

r2/„\ _ 71/T2 , tri/2 t;9 f dp 1 



M 2 (e)=M 2 + h5M 2 -h- ^- - „ (5) 

F(\ S?+ n f dp ■ / 2 ■ M2 ■ M 2 (e)-M 2 \ gh 2 ( f dp 1 \ 2 

S{e) = 6S+ 2 i< 1A 27 [y p +M + UWTW) + ^2{J y.wtw) (6) 

where M 2 , 5M 2 , SS and 8 remain finite when the cut-off is removed. The ambiguity in the 
choice of counterterms represented by SM 2 and 5S is resolved, as usual, by renormalisation 
conditions. We shall soon see that there is a natural way to do this in the present context. 
If F € [ip] is evaluated for a tp whose Fourier transform is non-zero only for momenta less 
than m it will reduce to a sum of local functionals of (p, 

F e [p] = JdxY: f J0 .^Mx) J °<P'(x) J1 --<P H (x) Jn (7) 



The expansion functions, (fj°(p'i 1 (p" : > 2 .., are related by partial integration so we can specify 
a linearly independent basis by insisting that the power of the highest derivative be at 
least two, and we will assume parity and (p —>■ —if invariance, restricting both the total 
number of f and the total number of derivatives in the expansion functions appearing 
in (|1|) and (|7|) to be even. It is important to note that (0) is not the same expression 
that would be obtained by acting with A e on the local expansion ([]]), because the former 
correctly includes differentiation with respect to the Fourier modes of if with momenta in 
the range rri^ < p 2 < 1/e absent from the second. Now if we scale the cut-off ( A ae W) [(p s ] 
extends to an analytic function in the complex s-plane with singularities only on the 
negative real axis, |3(J, the same is true of M 2 (se) and £(se), and consequently of the 



coefficents of the linearly independent expansion functions, fj ..j n in (^), so the contour 
integral 

WW = ^-l , ^e Xs V^sf J0 ., n (se) (8) 

Am J|s|=oo s 

can be calculated by collapsing the contour to a small circle about the origin and a contour 
along the cut on the negative real axis. The contribution from the circle about the origin 
is controlled by the small e behaviour of /j .. Jn (e)- As e — > this vanishes due to the 
Schrodinger equation, and in perturbation theory the Feynman diagram expansion gives 
an asymptotic expansion of fj ..j n (e) in positive powers of y/e. The inclusion of yn\s in 
(|8[) ensures that the contribution from the origin will be of order 1/A rather than 1/vA- 
For large \s\ the scaled field tp s is slowly varying and the scaled cut-off l/(se) is less than 
m so (A se W) [ip s ] can now be calculated by acting with A se directly on the local expansion 
of W, (U). Furthermore as the real part of A tends to infinity the contribution from the 
cut tends to zero due to the exp(As) factor, (the contribution from that part of the cut 
for which \s\ is large is again given by the local expansion and seen to be suppressed as 
A — > oo). Thus the Schrodinger equation leads to an infinite set of algebraic equations 
lim A _^ oo i i0J1 .. jn (A) = where 



T C(\\ *^^ ( U i -^0,2 A -Bo,o,2A 
I = -£{X) - n—= B 2 -\ j: h 



2 



7T V 3 10 J 

^ \ 3 10 



and 



u = % - sb 2Bi - n^ (ib* + »£ + °*e» + ..) 

I«,2 = \~ 1B 2 B (K2 - h^ (b 2i2 + 2B„, 4 A + i^^ + \ (9) 



5(A) = — / -e Xs VmVsS(s) = E h n S(\r (10) 

llll J|s|=oo S 

M 2 (X) = —[ —e Xs Vm\sM 2 (s)=M 2 + hM 2 (\) n (11) 

2ni 7| s |=oo s 



As the product se now plays the role of cut-off, rather than e alone, we have taken e 
to be finite and equal to unity. We will now choose renormalisation conditions. Note 
that the counter-terms only enter I and I 2 . If these are fixed then the above equations 
determine the coefficents Bj j n and the energy eigenvalue, £ , which are themselves finite 
as the cut-off is removed. Alternatively we could instead choose the values of two of these 
quantities, B 2 and £ for example, and then think of the equations 1$ = and I 2 = 
as determining the counter-terms. So we will take B 2 = —M/2, which is its classical 
value, and £ = as our renormalisation conditions. The advantage of imposing the 
renormalisation conditions on £ and B 2 is that we are free to solve © for the remaining 
Bj ,~,jn without first computing the A-dependence of the counter-terms which in a more 
general context can only be done in perturbation theory. 

The equations (|]) may be solved in the usual semi-classical approach in which we 
expand the coefficents as B = J2 h n B , by first ignoring the terms proportional to h. 
Although the resulting equations are quadratic in the Bj ^ n they are readily solved by 
starting with the coefficents of local functions of the lowest dimension and number of tp, 
giving at tree-level 

rjr f J I 2 1 12 1 *■ 112 *■ 1112 1 " IIII2 *■ 4 1 1 2/2 

W^ = Jdx[--<p --<p + -<p --<p +— V -- g(p +- g(p(p 



128 T r 256 T 1024" r r 256 r r 1024^ 

' 2 ////2 1 ^1 II 1112 1 '^ 12 1/12 "" "4 1 \ / 10 \ 



2048" r r 1024" rr r 2048" r r 4096" 



where we have chosen our mass-scale so that M = 1. Particular tree- level coefficents that 
will be of use are 

1(1/2" 

2 \ n 



Bi,,,^2 = -~ ( r ) (13) 



and the coefficients B\, B\ 2) B\ Q2) . . . , B\ ^ n=2 . To simplify the following formulae 
we re-name some of the coefficents. Firstly let B\ = bo/2, Bq o,j n =2 — b n ,n = 1,2,.. 
and B\ = c /6, B\ 2 , B\ Q2 , . . . , B\ Q Q , =2 = c n then the tree-level contribution to the 
equations 7 4 = 0, J 2 ,2 = 0, h,o,2 — 0, . . . , i2,o...o,2 = can be written as 

b d + 61 c = 

boC2 + Wei + b 2 c = 
b c 3 + b x c 2 + b 2 ci + b 3 c = 

etc (14) 

which, in turn, may be expressed as the vanishing of each coeflicent of z in 




IE 6 "*") (E C ^ m )- & 0C0 = 0, (15) 

We can solve for the c n in (|T5| ) to give 



R 1 

- D 2,0,.., 0,j n =2 



9_ 

16 



Rl Rl 

• D 0,2 ^O.O^ 

_1 5 d.2 
-1 



Rl Rl 

- D 0,0 v -Jn-i=2 - D 0,0,..J n =2 

Rl Rl 

- D 0,0,..,j„_2=2 - D 0,0,..,j„_i=2 

Rl Rl 

- D 0,0,--Jn-3 = 2 - D 0,0 v .Jn-2 = 2 



^0,2 



(16) 



We can find in a similar fashion the coefficients B\ = fo/15, B\ Q ^..g = f n ,n = 1... 
They are determined by the tree level equations J 6 = 0, / 4i2 = 0, -£1,0,2 = 0, etc. Redefining 
£>2 = \b , B\ = |co allows us to write those equations as 



Ev" +2 E& 




£/i*' - (co) 2 - 26 /o = 



(17) 



1=0 



where each coefficent of z must separately vanish. If we set (J2m=o b m z m )~ = Em=o PmZ m 
as well as (X)£L c nZ n ) = Y^=o lnZ n we can formally write, after substituting back the 
original values of b , c and / , 



1 1 n 

fn = ~ ((Co) 2 + 2b f ) /?„ - - E fain 
1 l fc=0 



1 



2 V 27648 



Pn + E Pkln- 



k=0 



n > 1 . Using the formulae for inversion and product of power series from the mathemat- 
ical literature |32| , we can calculate all the B\ Q Q ^ =2 



The order-fr corrections are obtained by substituting the tree-level results into the 
previously ignored order-fo term in the Schrodinger equation and treating this as a per- 
turbation to the classical equation. We want to use this to show that our large-distance 
expansion correctly gives the short-distance behaviour as contained in the divergent mass 
and energy subtractions S(X) and M 2 (X), which occur only in J and I<i- We first study 
I2. Using ( 16|) we get the 0(h) expression 



7»(A) 



^t + 2Bi 



V\ ( 1_ _A_ 



A 2 



9 



A V 16 



192 + 1280 



5A 3 



43008 



+ .. 



(19) 



This vanishes when A — ► 00, but we get a good approximation if we truncate the series 
and take A as large as the truncation will allow, i.e. small enough for the first neglected 
term to be insignificant. Since I(X) is of order 1/A for large A the accuracy of this 
approximation is greatly improved if we perform a further contour integration, amounting 
to a re-summation of the series in A. Observe that substituting A = 1/yfs in I(X) gives a 
function that is analytic in s with a cut on the negative real axis that we wish to evaluate 
as s tends to zero from real positive values, so we define 



j( A ) = J_ f ^e A2s v / ^7/( S - 1 / 2 ) 

2m J\ s \=oo s 

for which lim^oo I(X) = 0. Thus 



m = 6 -^P^ + 2B^ 9 s(x) 



(20) 



(21) 



where 



5(A) 



^ 1 161X3/4) 



A v / 2r(3/4) 



96tt 



+ 



A 2 



960r(3/4) 



A 3 x/2r(3/4) 
5376tt 



(221 



The terms in 5(A) now decrease more rapidly than the corresponding terms in /(A). 
Since I{\j\fs) behaves asymptotically as \fs for small s this re-summation has the effect 
of eliminating the leading term so that /(A) is now of order 1/A 2 . Further re-summations 
are only efficacious given a sufficent number of terms in the truncated series for the extra 
gamma-functions in the coefficents to be noticeable. 
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Figure 1: The mass subtraction 

In fig. flip we plot, the series S(X) truncated to 13 terms, —M 2 (X) n /(2g), their sum, 
and the limit of this sum as A — > oo, (which we obtain exactly in the next section as 
-1/(8tt) ~ -0.0398). Clearly neither 5(A) nor -M 2 (X) h /(2g) are constant for large 
A but their sum is, to a good approximation for A > 2. This shows that our large- 
distance expansion correctly reproduces the short- distance effects encoded in M 2 (e) h . The 
departure from this constant value for A > 11 is due to the error involved in truncating 
5(A) to 13 terms. If we denote by S n the series truncated to n terms minus M 2 (X) h / (2g) 
then in fig. (|P we have shown S n for n = 4, 5, 8, 9, 12, 13. 

Each truncation provides a good approximation to 5(A) up to a value of A which is 
large enough for the highest order term to be a significant fraction of the whole. Taking 
this to be one per cent gives an estimate of 5(oo) with an error that ranges from three 
per cent (five terms) to half a per cent (13 terms). 
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Figure 2: Truncating 5(A) 



To check that our large distance expansion correctly reproduces the energy subtraction 
we need the 0(h) part of W[ip] that is quadratic in ip. We obtain this from the equations 
Io,2 = 0, Jo,o,2 = 0,.., having imposed the renormalisation condition B 2 = —M/2. We 
use the re-summation described earlier, truncate the series in A so that they include 
contributions from coemcents of functionals of ip of dimension less than 26, and take A so 
that the last incuded term is one per cent of the value of the truncated series. We also 
use Stieltje's trick of halving the contribution of the last included term to improve the 



accuracy of the approximation |33[]. This gives the estimate 



W. 



n 9 



,(5)2 



I dx 1 6.6 V 2 - 6.02y/ /2 + 5.40V" 2 - 4.9V" 2 + 4.54y? ( 
-4.24<^ 2 + 4.0M 7 ) 2 - 3.79</ 8 > 2 + 3.58<^ 2 - 3.34<^ 10 ) 2 + .) (23) 



In the next section we obtain W 7 ^ exactly. Rounding the exact results to three signif- 
icant figures gives 



W. 



h _ 9 



,(5)2 



. ( I dx ( 6.63<// 2 - 5.97y/ /2 + 5.33y/" 2 - 4.84y/'" 2 + 4.45y? ( 
-A.Ucp^ 2 + 3.89^ 7 > 2 - 3.68<^ 2 + 3.50<^ 9 > 2 - 3.34<^ 10 ) 2 + .) (24) 



which shows that our approximate results are good to a few per cent. 



Energy Subtraction 
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Figure 3: The energy subtraction 

Figure @ shows the effect of substituting this estimate into the 0(h 2 ) contribution 
to Jo- The top curve, A, is the estimate of the re-summation of the series in A, whilst the 
bottom curve, B, is the 0(h ) contribution to the re-summation of £(A) evaluated using 
(^D with 5M 2 = g/(4:7r). Neither of these curves appears to tend to a constant for large A 
whereas their sum, represented by the middle curve, C, provides a good approximation to 
a constant value for A larger than four until A is sufficently large that the approximation 
of the infinite series by just ten terms breaks down. The straight line in the figure is the 
value 0.0052 which would be obtained by truncating the series at fifty terms using the 
expression for Wf we find in the next section. 

Having seen that our large distance expansion successfully reproduces the short- 
distance effects contained in the counter-terms of the Hamiltonian we turn to the one-loop 
evaluation of the Bj ^j n coefficents corresponding to higher numbers of fields. Begin with 
the coefficents of local functions containing four fields. Rounding the tree-level result to 



,M 



three significant figures gives 

W\ = j^ J dx( -10. V + 15.6^V 2 + 3.9 V 4 
-7.81</?V 2 - 30.3^'V" 2 - 11- 7w" 3 + 4.88^V'" 2 - 22. 7y/' 4 

+36 W" 2 + 40.0 W 'V'" 2 - 3.42^ V"" 2 + •■) (25) 

Estimating the 0(h) contribution in the same way that we estimated W£ gives 

W? = j^ J dx U.02<p* - 20.(ty V - 7.96/ 
17.W" 2 + 83.8^'V" 2 + 37.6 W " 3 - 15.6<p 2 (p m + 87.7y/' 4 

-12 W" 2 + -164^^'V'" 2 + 14.0^V"" 2 + ••) (26) 

There are two things to note. Firstly there is a proliferation of local functionals of the 
same dimension and number of <p as these increase. So, for example, there is a unique 
local functional with just two (p for any dimension, but there are two hundred and seven 
with twelve ip and dimension twelve. Secondly the ratio of the 0(h) corrections to any 
two coemcents of functionals containing the same number of ip and the same dimension 
is approximately the same as the ratio of the tree-level values. For example the ratio 
of the 0(h) coemcents of <p> 2 <p>" 2 and cp' 4 is —17.4/7.96 ~ —2.19.. whereas the ratio of 
the corresponding tree- level values is exactly —2. Given that our estimate is probably 
only good to a few per cent it is not clear at this stage whether the one-loop ratios are 
exactly equal to the tree-level ratios, but we will investigate this with greater accuracy in 
the next section. We will now compare these results with those obtained by solving the 
Schrodinger equation without first expanding in terms of local functions. 

3 Direct Semi-Classical Solution 

It is straightforward to solve the Schrodinger equation lim e j F e [<p] = without resorting 
to the local expansion, at least for low orders of an expansion in ip and h,^. This turns 
out to be remarkably efficent compared to the Feynman diagram expansion which we 
describe in the next section. Expand W[(p] as 

oo . 
W [v] = Yl J dp 1 ..dp 2n <p(Pl)--<p(P2n)^2n(Pl, ■■,P2n)$(Pl + -+P2n) (27) 



n=l 



where the T are unknown functions. Then we can write A £ W[<£>] = J2 Ar 2n where Ar 



2/i 



IS 



/ 2irdq / dp 3 ..dp 2n 2n(2n - l)(p(p 3 )..<p(p 2n )T2n(q, -q,P3~,P2n)$(j>3 + ■■ + P2n) (28) 

Jq 2 <l/e J 



and 



/^(f) 2 = gr 2 „or 2 ,„ (29) 



where T 2n o r 2m is 



8nmVT / dp 2 ..dp 2n dk 2 ..dk 2m ip(p 2 )..lf(p 2n )ip(k 2 )..ip(k 2m )T 2n (-(p 2 + •• +P2n),P2,--,P2n) 

xT 2 rn(-(k 2 + .. + k 2m ), k 2 , .., k 2m ) 5{p 2 + .. + p 2n + k 2 + ..k 2m ) (30) 

Expanding the T in powers of h as T 2n = J2 ^ m ^2n anc ^ ignoring oidei-h terms in the 
Schrodinger equation gives the tree-level result 

r* o Tl + 2Tl o T\ + 2Tl o Tl + T\ o T\ + .. = J dx (u/ 2 + M 2 <p 2 + ^ A \ (31) 

The term quadratic in if is 

TloTl = J dp8n(rl(p,-p)) 2 v(pM-p) = f ^(P 2 + M 2 )p(p)p(-p) (32) 

so if we take the negative root for normalisability of the vacuum functional we get 

! y/p? + M* _ W(p) 

1 o — . — 

Using 

TloT 2n = - I dpi..dp 2 n<fi(Pl)- t f{P2n) {j2^{Pi)j ^2n{Pl, ■; P2n)5{Pl + •■ + P2n) (34) 

r£ori = !|dV. ( 35 ) 

^ives, for pi + .. + p 4 = 0, 0, 

rl(pi,--,p 4 ) = ~ / w ,,w / x , — ; — /— \V ( 36 ) 

(27ry{4\)(uj{p 1 ) + .. + uj{p 4 )) 



2 - A = — 1-- ( 33 ) 

47T 47T 



111 



The terms of higher order in (p, for which there are no contributions from the potential 

r 1 or 1 

n+m=const x 2n 2m 



give J2 n +m=const ^ln ° ^Im = which can be solved recursively as 



ri.(pi, ..par) = 

47T r_1 r 

^2r / N E n ( r + l ~ n ) S { r 2n(-(P2 + •• +P2n),P2, -,P2n) 

xr 2(r+l-n)(-(P2n+2 + •• + J?2(r+l-n)), P2n+2, ..,£>2(r+l-n)) } (37) 

where S symmetrises the momenta. Expanding T\ and T\ in positive powers of the 
momenta reproduces (|i~2| ) as it should since no re-summation is involved in either approach 
to the tree-level result. 

The order-fr contribution to the Schrodinger equation is 

E rL o r|L + E Ar L + [** (2£ h - (m 2 ) V 2 ) = o. (38) 

n,m n 

10 



The term quadratic in tp gives the limit as e — > of 

5M 2 + M 2 (e) h , x „ fi/ N o rWi do , N 

- 2w(p)r*(p, -p) + -M- -7-^-7-r = (39) 



4tt " v ^ y ^' ^ 32tt 2 J-i/y/eu(q)+u(p) 

which can be solved for p ^ as, |l|], 

q f°° f 1 \ <5M 2 

r5(P,-P) = ^/ o ^ UmM^) + ^p)) J"^w (40) 

- » sinh -! fiL) _ J£L (41) 



32ir 2 p \MJ 8iruj(p) 

and for p = we get r| = g/(32ir 2 M) —SM 2 /(8irM). The renormalisation condition that 
fixes Bi at its classical value requires that r|(0, 0) = 0, which determines SM 2 = g/(47r). 
Setting p = in (|39|) and taking the limit e — > is meant to yield the same as taking the 
limit A — > oo of /q(A) i n dUD when we identify r 2 (0,0) = S 2 /(27r). This gives the value 
—g/(8ir) quoted earlier that agrees well with the large A behaviour of 5(A) — M 2 (X) h /2. 
More particularly S(X) should be obtained from the large e expansion of 

m = — ( r * q t , (42) 

167T Jq 2 <l/e (jj(q) + Uj(p) 



by applying the two contour integral re-summations, giving 

A 

'|s|=oo ^\s\- 



-—I dss~V 4 e x2s [ dss-^e'/^His) (43) 

A-K J\s\=oc J\s\=oo 



which does in fact coincide with fl2"2"|). Since the large A behaviour corresponds to small e 
we can use (O) to investigate this. Thus for small e 



dq[ . . . - . — = v / M- 1 + e- v ^^v / M-v / ^+^p= ^=a+» ( 44 ) 



which leads to the power law corrections to the large-A behaviour described earlier. Ex- 
panding (|4i"l) in positive powers of p 2 leads to the exact results for the W 2 quoted earlier 

J 27^^' ~ p ^(P^(~ p } = 
g f , /V 2 3y/' 2 15/" 2 35</" 2 315y?( 5 ) 2 
7i J \48 160 896 2304 22528 

693<^ 2 lOOl^ 2 6435^® 2 109395^ 9 ) 2 230945<^ 10 ) 2 \ 
+ 53248 + 81920 + 557056 + 9961472 + 22020096 + "J ^ ' 

The 0(h) contribution to the part of W[(p] that is quartic in tp is obtained from 

2T\ oT 4 s + 2T\ o Y h 2 + AT\ = (46) 
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which leads to 

r4(Pl,~.P4) = .,,_,,„* 4 1 N S 



^ s f f 00 dq ( 1 

(27r) 3 4!7rEi^fe) [Jo 2co(q) + E^(Pi)\ oo(q)(u(q) + u( Pl )) 



3 

4 



u;(g) +w(pi) +uj(p 2 ) +u(q + p 1 + p 2 ))(uj(q) + uj(p 3 ) 4-^(p 4 ) + io{-q + p 3 + p 4 )) 

1 1 



2w(Pi)Eiw(pi 



(47) 



Expanding in positive powers of the momenta p 1; ..,p 4 and integrating over g numerically, 
(using MAPLE), leads to 

W * = mOO J dx ( 3 - 973 ^ ~ 19 - 45 ^V 2 - 7.961<p /4 

16.27(^V" 2 + 85.78^'V' 2 + 33.66^" 3 - 14.15^V" 2 + •• ) • (48) 

From this it is clear that our previous estimate was quite good, but that the observation 
that the ratios p 1n ,- = B\ ■ /(qB° n ■ ) are the same for coefficents of functionals of 
the same dimension and number of fields is only approximate, since 

p 4 = -0.03814 (49) 

p 2 , 2 = -0.1245 (50) 

p ,4 = -0.2038, p 2 ,o,2 = -0.2082 (51) 

Po,2, 2 = -0.2834, p 1A3 = -0.2872, p 2 , ,o, 2 = -0.2898. (52) 

These ratios may be explained by observing that the dominant contribution to ( fJ7|) comes 
from the last term in the braces, which itself originates in the mass renormalisation. hSM 2 
was fixed by imposing B\ = 0, but if instead we had taken hSM 2 = then this term 
would have been absent. The effect of this choice on the ratios can be calculated using 
dimensional analysis and gives 

p 4 = -0.00165 (53) 

p 22 = -0.00513 (54) 

p 0]4 = -0.00486, p 2 ,o, 2 = -0.00926 (55) 

Po,2,2 = -0.00488, pi, ,3 = -0.00868, p 2 , ,o,2 = -0.0113. (56) 

The advantage of this choice is that the one-loop corrections to the coefficents Bj 0i „j n 
for functionals containing four fields are now significantly smaller. The same is true for 
the coefficents corresponding to two fields, with the exception of B 2 . This suggests that 
a more effective choice of renormalisation condition which would reduce the size of the 
one-loop corrections, would be to fix i? 4 at its classical value, rather than B 2 . 
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4 Sinh-Gordon Model 

From the standpoint of perturbation theory <p 4 -theory is 'close' to a theory that is quite 
special, namely the Sinh-Gordon theory which has an infinite number of conserved quanti- 
ties that imply the absence of particle production, it is interesting to calculate the vacuum 
functional for this case. The potential of the Sinh-Gordon theory may be taken to be pi 



M 2 + MM 2 ( (3 2 r 

V = - 2 cosh^exp^--^ 



l/Ve dp 



Lj(p) 



r2 , fcJf waw 1 , <? , 0V , 0V , \ A P 2 f 1 ^ dp 



Apart from the replacement g — > M 2 (3 2 the Sinh-Gordon potential leads to the same 
expressions for the tree- level values of r^T^ and the one-loop result Y\. The tree-level 
Tg is modified by the <p 6 term in the potential 

/3 4 M 2 
r ^ ^ " 6!(2^EXS (58) 

this, together with the <5M 2 <p 4 term in V modifies the one-loop value T^ 



ri-r}- ,-^.^.,., / "» „.,.,. U.,., -^ +^ <») 



/5 4 m 2 ( r a ( l l \ 1 

so that for the sinh-Gordon model 

rrrn P A f , I <P* 5tt - 22 2 /2 2275tt - 8952 /4 651tt - 2768 2 //2 
4 7T i ^384 1280 Y ^ 860160 r 172032 Y Y 

1041 705vr - 4243072 /2 //2 689535tt - 2920448 ,„ 



41287680 r r 82575360 

13905tt - 58624 2 ///2 N 

3932160 ^ ^ 



/3 



i 



/ ds(8.2893^ 4 - 15.647<pV 2 - 6.6791<p /4 



10000 

+13.3743^V' 2 + 74.818<pV' 2 + 29.0731<p(p" 3 - 12.0941^ 2 ^'" 2 + .. J . (60) 

The ratios of the one-loop coefficents to their tree-level values for this model are 

p A = -0.07958 (61) 

p 22 = -0.1001 (62) 

p ,4 = -0.1710, p 2 , ,2 = -0.1712 (63) 

Po,2, 2 = -0.2471, p 1A3 = -0.2481, p 2 , , ,2 = -0.2477. (64) 

Note that again the ratios are approximately the same for coefficents of functionals of the 
same number of fields and dimension, however this cannot be explained away as simply 
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the effect of mass or coupling renormalisation. If we had taken 5M 2 = we would have 
obtained 

p 4 = -0.1194 (65) 

p 22 = -0.06035 (66) 

p 04 = -0.05162, p 2A 2 = -0.05183 (67) 

Po,2, 2 = -0.04820, ,01,0,3 = -0.04915, p 2 ,o, 2 = -0.04875. (68) 

which are approximately constant for coefficents of functionals of the same number of fields 
and dimension. These coefficents would, with a change of sign, apply to the Sine-Gordon 
model as well. 

5 Feynman Diagram Expansion 

We will now describe how the results of the previous section are obtained within the more 
conventional approach to field theory based on Feynman diagrams. There are several 
ways to represent the vacuum functional as a functional integral, the most convenient for 
our purposes is the representation, 0, 

y[(p] = f V4>e- SE[(l,]+ -f ' dx ** (69) 



where <ft vanishes on the surface t = 0, which is the boundary of the Euclidean space-time 
t < in which the field lives. Therefore, W[ip] is a sum of connected Euclidean Feynman 
Diagrams in which ip is the source for on this boundary, the only major difference from 
the usual Feynman diagrams encountered in free space is that the propagator vanishes 
when either of its arguments lies on the boundary. Such a propagator can be obtained 
using the method of images as 

G(x,y) = G (x,y) - G (x,y) = G„(x,y) - G„(x,y) (70) 

where x = (t,x), x = (— t, x), similarly for y and y, and Go( x 5 y) is the free space 
propagator. Developing the Feynman diagram expansion yields the tree level diagrams 
up to <p 6 which are shown in figure (HI). 





(a) (b) (c) 

Figure 4: Tree Diagrams up to <p 6 
The Feynman rules in the coordinate space for the diagram (b), for instance, yield 

gs f J d Xl ..dx A d vvMMx*) — Wi Wi Ws ^"— (71) 
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where Sf is the symmetry factor associated to the diagram, the times £1, ..,£4 are all set 
to zero, and the time-like component of v is integrated over negative values only. When 
a propagator ends on the boundary it is equal to twice the free space propagator, so that 
in momentum space this is proportional to 

/ dpi..dp4 dq 1 ..dq 4 0(p 1 )..<p(p 4 )6(p 1 + ..+pa) %i + .. + &) J[ 2 , 2 J , . , 2 (72) 

The momenta <7i, .., 94 are the time- like components which, when integrated over, give rise 
to 5- functions that impose the vanishing of £1, ..,£4 in the configuration space represen- 
tation (|7T|). In the momentum-space representation these integrals are readily performed 
as contour integrals leading to our previous expression (p5G|) for T\. Diagram (a) gives 

j^jdpdq <p(p)<p(-p) p2 + ^ + M2 (73) 

The integral over g is divergent. Formally we can write it as 

J ^ v>(p)#-p) (^(0) - vS + M 2 ) (74) 

The origin of 5(0) is in the construction of the original path integral representation of W 7 , 



30]] . Rather than the functional integral (p9| ), we should start with 

V[cp] = (<p\Q) = (D\e i f* v dx \0) (75) 

(where the state (D\ satisfies (D\ tp — 0) . As operators tt = cf) , but in the passage from 
([75|) to the functional integral % is represented by <fi plus terms coming from the time 
derivative acting on the T-ordering because the functional integral represents T-ordered 
products. This leads to terms like / dx(p 2 5(0), because this is local it may be cancelled by 
an equal but opposite counterterm, which amounts to simply discarding the divergence. 
Alternatively, and perhaps more satisfactorily, we can deal with this divergence by placing 
the source terms itip not at £ = but at small, distinct times £, and finally taking the 
limit £j — > . This leads to an extra factor of the form e tqe in ([73]) which regularizes 
the divergence and enables the integral to be done by closing the g-contour in the upper 
half-plane, yielding a finite result in limit as e — > 0, so that we end up with our earlier 



expression (|33D . 

When neither end of a propagator is on the boundary, £ = 0, the image charge breaks 
energy conservation leading to more complicated expressions. However for the tree-level 
diagram (c) these may be combined into an expression proportional to 

(yr dpjdqi <p(qj) pA S(pi + ..+p 6 ) 5(q 1 + .. + g 6 ) ( 

\\ l P* + q? + M 2 ) (p 1 +p 2 + p 3 y + (q 1 + q 2 + q 3 )'> + M2 [ ] 

which, after integration over the <& gives the same as fl3"7D- 

The diagrams represented in figure ([5]) enable us to calculate the 0(h) correction for 
the coefficients of 2 and 4 fields . 
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(hi 





(d) (e) 

Figure 5: Two and four field one loop diagrams 
From fig|)](a) above we can work out the first order correction to B 2 as 

s f g / dx 1 dx 2 d 2 z <p(xi)(p(x 2 ) — ^r — G(z,z) — -- 1 — (77) 

J Ot\ ot 2 

This is divergent due to the factor G(z, z). Since the propagator does not touch the 
boundary, it has an energy non-conserving contribution from the image represented by 
the exponential in 

c(77 , fdpdq (1-e 2 ^) 

where t is the time-like coordinate of z. We regulate the integral by restricting the 
space-like component of momentum, p 2 < 1/e, just as we did for the Laplacian in the 
Hamiltonian. The divergence is then cancelled by the counter-term represented by fig.|5](b), 
which is due to the order-fr terms in M 2 (e), (|j). Taken together these diagrams give T 2 
as in (|T]). 

The first order correction Tj is given by the sum of diagrams |5](c), |5|(<i) and |5](e). As 
before, diagram (e) is the counterterm diagram associated to the bubble appearing in (<f). 
In momentum space ^(c) gives 

2 (c) 

—^— / dkdp 1 ..dp i ip{p\)..<p{pi) 5(pi + .. +pi) 

Ei^(Pi)+^(fc)+^(0 

£i u>(p0(a;(p3) + u ^) + w (*0 + w (0)( w (Pi) + W (P2) + w(fc) + w(0) ' 

(E^fe) + 2^))(Ei^fe) + 2c(/))' (79) 

where Z = £>3 + p^ + fe. whereas, for the sum of the other two diagrams, (with SM 2 = 0), 
we have 

^ 'dkd Pl ..dp t „„ , f (Pl) ;,^ P ' ) / (Pl ^ 4 + f 4 ' jrrr (80) 



t 4 / "o>(A)(w(j)4)+wW)EtwW(E;wfe)+2a>(«:)) 

The symmetry factors Sf and s\ are respectively equal to 1/16 and 1/12. If we expand 
in powers of the momenta pi, --,P4 the loop integration over k can be done to reproduce 
the results for W± given by (f£8p. 

The extra (p 6 term in the potential of the Sinh-Gordon model generates the diagram 
fig. (|). Its analytic expression after subtracting the counterterm reads 

1 f j, J j £(pi)..£(p4)<S(pi + ..+P4) , 01 v 



/3 4 s f77T~)3 / dkdpi..dp 4 



{2-nfJ ^" ^ cv(k)(2u(k)+T,M Pi )) 
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where s/ for this diagram is 1/96. This gives the modifications to the cp A results described 
earlier. 





(a) (b) 

Figure 6: 6 point interaction for the Sinh-Gordon Model 

6 Reconstructing the Vacuum Functional 




Figure 7: Tree- level vacuum functional: T\(p, —p) 

We will now use the above results for the <p 2 part of W[ip] to illustrate how the vacuum 
functional can be reconstructed from its large distance expansion. The tree-level contri- 
bution has been treated in detail in fl30f . Applying formula (||) to / dp(p(p)ip(—p)T2(p, —p) 
amounts to the expression 



_ , . , If ds 
T 2 (P, -P) = hm — - / -. 

A^oo 2m J\s\=oo s — 1 



X(s-l) 



\ZsT 2 {p/Vs, -p/y/s) 



(82) 
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Since \s\ is large on the contour we can use the local expansion T 2 (p, —p) 
Shifting s we get 



2^0 ®"n.P 



2n 



lim 



J\s\=O0 S r. 



P 



2n 



(s + D 



= lim S(p, A) 

A^oo 



(83) 



Expanding the (s + 1) factors in powers of 1/s enables the integral to be done, yielding a 
power series in A. For example, at tree-level 



1 

2tu 



ds 



,*(«-!) 



|s|=oo S 



p 2 + s = J2 



(_)n+l A n-l/2( 1+p 2)n 

n!(2n-l)v^r 



M) 



This series converges for all positive A. We get an approximation by truncating the 
expansion by including terms up to and including A^ -1 / 2 , say. This requires a knowledge 
of the local expansion only up to terms in (ip( N '(x)) 2 . To demonstrate this approximation 
we have plotted in fig. (0) the series ( [Si]) truncated at iV = 12, Sw The value of A 
is chosen so that the last term included is one per cent of the value of the series. We 
have also plotted \Jp l + 1, and the expansion of \fp 2 + 1 in powers of p 2 , C truncated to 
fourteen terms. The full series fails to converge for p 2 > 1, and this is reflected in the 
fact that the truncated series ceases to be a good approximation for p 2 > 1. However, 
S12, which is a resummation of this series is a very good approximation for a much larger 
range of momenta. 
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Figure 8: One-loop vacuum functional: T%(p, —p) 

The one-loop correction, T§, may be treated in the same way. In fig. (^) we have 
plotted arcsinh(jo)/p. The small p expansion, C is again only good for p 2 < 1. Our 
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approximation that re-sums this series, Su provides a good approximation only over a 
slightly larger range. The accuracy of this approximation is greatly improved by further 
re-summations, just as we did for the Laplacian. Let us define the re-summation operator 
acting on a function of A and p to be 

R.S(X,p) = ^-[ ^e Xs S(s-V\p). (85) 

2m J|s|=oo s 

Then the curve R\2 shown in fig. (|j) results from applying R twice to S12, and provides 
a good approximation to arcsinh(p)/p for values of p up to about p — 5. Since the 
effect of applying RP to a term in S12 that is proportional to A n is simply to divide it by 
T(n/2 + l)T(n/4 + l)..r(n/2 p + 1) further applications would have no significant effect 
when we take just twelve terms in the expansion. 

7 Conclusions 

The purpose of this paper has been to test an approach to quantum field theory in 
which states are constructed in the Schrodinger representation from their large-distance 
behaviour. The vacuum functional is expanded as a series of local functionals. Truncat- 
ing this series at some convenient order leads to an approximation scheme. For scalar 
if 4 theory in 1+1-dimensions we compared the results obtained by solving the form of 
the Schrodinger equation that applies in this approach with those obtained from a semi- 
classical expansion, and from the Feynman diagram expansion of the wave-functional. We 
have found that the expansion coeficents agree to within a few per cent when we truncate 
these series at about ten terms. Also we found that the known counterterms that contain 
information about short-distance effects are correctly reproduced by our large-distance 
expansion to a similar accuracy. 

We also found a curious simplification that occurs in the vacuum functional of the 
Sinh-Gordon and Sine-Gordon models. The ratios of coefficents of the one- loop corrections 
to the coefficents of local functionals containing four fields to their tree-level values are 
approximately the same for functionals of the same dimension. 
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